一类Marangoni对流边界层方程的近似解析解

利用Adomain解析拆分和Padé逼近方法对由Marangoni对流诱发的层流边界层问题进行了研究, 提供了一种求解边界层方程的解析分析方法. 得到了问题的近似解析解并对相应的流动及传热特性进行了探讨. 本文所提出的思想方法可以用于解决其他科学和工程技术问题. 关键词： Marangoni对流 非线性 Adomain拆分法 近似解析解     

A dynamic multiscale lifting computation method using Daubechies wavelet

An important property of wavelet multiresolution analysis is the capability to represent functions in a dynamic multiscale manner, so the solution in the wavelet domain enables a hierarchical approximation to the exact solution. The typical problem that arises when using Daubechies wavelets in numerical analysis, especially in finite element analysis, is how to calculate the connection coefficients, an integral of products of wavelet scaling functions or derivative operators associated with these. The method to calculate multiscale connection coefficients for stiffness matrices and load vectors is presented for the first time. And the algorithm of multiscale lifting computation is developed. The numerical examples are given to verify the effectiveness of such a method.     

Mixing mechanisms of vortex ring formed by gravity slumping motion

 This paper is concerned with an experimental investigation of the mixing inside the vortex ring formed by the gravity slumping motion of a dense cloud in a less dense atmosphere. The dynamics of the spreading and instantaneous structures of the turbulent flow were examined by visualization, single and multi-point measurements of velocity and concentration for two heavy gases, carbondioxide (CO₂) and dichlorodifluoromethane (CCl₂F₂), in a configuration in that heavy gas, initially trapped in a reservoir, was released with the rise of a shutter into calm air of a sector-shaped dispersion channel. Visualization of the cloud as a whole showed a spreading motion in which an advancing frontal structure was followed by a stratified flow with a layer of dense fluid of higher velocities near the wall and, on top of it, a layer of dilute fluid whose concentration is controlled by the mixing mechanisms within the head. During the course of spreading, there was always a phase in which the head attained to a constant speed of advance, which occurred as 0.13 m/s for CO₂ and 0.48 m/s for CCl₂F₂. It was interesting to observe for CO₂ that the phase of constant speed took place in between two acceleration phases; the former was due to the initial slumping of the cloud at the exit of the reservoir, and the latter was attributed to the collapse of the head on the transition to the passive dispersion phase. Instantaneous two-dimensional velocity field, measured with particle image velocimeter (PIV), showed that the cloud overran the ambient air which caused the approaching dense fluid deflected away from the wall with significant vertical velocities and downstream-moving separation, and the air trapped under the head resulted in the density inversion which introduced further intricacy to the turbulent structure of the head. Instabilities at the upper free shear layer due to density and velocity discontinuity rolled into periodic array of vortices which engulfed a considerable amount of air as they were convected backwards over the head, but the incorporation of heavy and light fluids was completed with the appearance of microscales after the collapse on the stratified layer. Analyses of the cloud head at different downstream locations also revealed that its size remained unchanged when the speed of advance was constant, allowing the rate of change of the cloud volume being modeled with the rate of spreading. Contours of concentration obtained from digitized PIV pictures confirmed the kinematic features of the mixing revealed by the velocity field and that the concentration values within the large structures were higher than those at the upper part of the stratified layer. Motivated by the experimental observations, a semi-empirical analysis was presented to describe the results and based on local values of the Richardson and Reynolds numbers. Received: 4 October 1995 / Accepted: 4 July 1996     

Oblique reflection of shock waves from rigid porous materials

A study to determine the general gas dynamic behaviour associated with the impact of a shock wave on a porous wedge has been undertaken. A number of interesting features are noted. The pattern of wave reflection is shown to be significantly affected by the inflow of gas into the wedge. This has the effect of reducing the triple point trajectory angle for cases of Mach reflection and for strongly reducing the reflection angle in regular reflection. The permeability of the wedge has a significant effect on the strength of the reflected wave and in some cases this wave can be attenuated to the extent that it is almost eradicated. Pressure measurements taken under the wedge are characterized by oscillations which are of similar shape, for a given wedge, over a range of shock wave Mach numbers. It is shown that the wave transmitted into the wedge is attenuated to varying degrees depending on the material properties, and that for weak incident waves the mean propagation velocity can be less than the sound speed in the pore fluid. Photographs taken using a specially constructed wedge which allows the transmitted wave to be visualised, show that the transmitted wave is nearly plane.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Electrically forced microthreading of highly viscous dielectric liquids

The behaviour of three high viscosity (4875, 12 125 and 58 560 mPa s), dielectric liquids was investigated at flow rates of 10⁻¹⁰, 10⁻¹² and 10⁻¹⁴ m³ s⁻¹ and the applied voltage range 6–15 kV. In these experiments, due to the low electrical conductivity of the liquids (10⁻¹³ S m⁻¹) and therefore the ensuing high electrical relaxation time, classical electrohydrodynamic atomization conditions are not satisfied. Only dripping and unstable jetting were observed at 4875 mPa s. A transition from no jetting to stable microthreading was observed for the 12 125 and 58 560 mPa s samples. The relics accompanying the transition were found to change from discrete droplets to a continuous filament. Stable microthreading, which generates uniform filaments, was obtained for the 12 125 mPa s sample at flow rates 10⁻¹⁰ and 10⁻¹² m³ s⁻¹ and in the case of the 58 560 mPa s sample at all the flow rates investigated. The high viscosity assisted stable microthreading with the filament diameter decreasing with increasing applied voltage and more dramatically decreasing with reducing flow rate.     

The effect of viscosity on the free stream in transonic flow over a corner point

An asymptotic analysis of the system of Navier–Stokes equations for describing the flow which arises from the subsonic free stream in the neighbourhood of the vertex of a convex corner with curvilinear generatrices is presented for Reynolds numbers approaching infinity. It is assumed that, in limiting non-viscous flow, the subsonic free stream reaches the velocity of sound at the vertex of the corner and, in the first approximation, is described by the Vaglio–Laurin solution. It is shown that the flow can have a different form depending on the value of the pressure gradient, which is formed in the neighbourhood of the corner point. However, irrespective of the steady form of the flow, as a result of the interaction of the Vaglio–Laurin flow with the boundary layer, the latter induces perturbations in the outer flow, which “rounds off” the vertex of the corner when there is a transonic flow around it.     

Global and local optimization using radial basis function response surface models

The focus of this paper is the optimization of complex multi-parameter systems. We consider systems in which the objective function is not known explicitly, and can only be evaluated through computationally intensive numerical simulation or through costly physical experiments. The objective function may also contain many local extrema which may be of interest. Given objective function values at a scattered set of parameter values, we develop a response surface model that can dramatically reduce the required computation time for parameter optimization runs. The response surface model is developed using radial basis functions, producing a model whose objective function values match those of the original system at all sampled data points. Interpolation to any other point is easily accomplished and generates a model which represents the system over the entire parameter space. This paper presents the details of the use of radial basis functions to transform scattered data points, obtained from a complex continuum mechanics simulation of explosive materials, into a response surface model of a function over the given parameter space. Response surface methodology and radial basis functions are discussed in general and are applied to a global optimization problem for an explosive oil well penetrator.